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Abstract 

At first, we generally investigate the short-time decoherence of a qubit 
nonlinearly coupling with a bath. The measure of the decoherence is 
chosen as the maximum norm of the deviation density operator. Then 
we concretely investigate the Josephson charge qubit (JCQ) model. It is 
shown that at the temperature T ~ SOmA", the loss of fidelity (due to 
decoherence) of the JCQ is bigger than the DiVincenzo low decoherence 
criterion. The decoherence will decrease with the decrease of the experi- 
mental temperature. When the temperature decreases to T ~ Q.?,mK the 
DiVincenzo low decoherence criterion can be satisfied. 

PACS numbers: 03.65.Ta, 03.65.Yz, 85.25.Cp 

Keywords: Short-time decoherence, Josephson junction, quantum com- 
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1 Introduction 

David DiVincenzo put forward a low decoherence criterion for the candidates of 
quantum computing hardware to be satisfied W. An approximate benchmark 
of the criterion is a fidehty loss no more than ~ 10~^ per elementary quantum 
gate operation. In general, the loss of the fidelity results from quantum leakage 
as well as decoherence. The Josephson junction is considered to be a promising 
physical realization of the qubit. It has been shown that the quantum leakage 
of Josephson charge qubit (JCQ) is not serious So it is very interesting to 
investigate the decoherence of the JCQ. To perform the quantitative study of 
the decoherence for a qubit, in general one needs to solve the quantum dynami- 
cal problem of the qubit coupling with its environment. However, the solutions 
of the coupling system can not be obtained easily, usually a kind of approxima- 
tion must be appealed. The Markovian approximation j2] is the most familiar 
one. This approximation has been used to study the decoherence of qubits in 
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last years [S]. But it is not a suitable tool for researching the short-time 
decoherence of the qubits [S]. However, the short-time processes is very impor- 
tant because almost all of the gate operations of the quantum devices are short 
time. Fortunately, a short-time approximation scheme for the reduced density 
matrix has been developed recently 0. Another core of the issue to investigate 
the decoherence of qubit is to make certain what the model of the environment 
and the coupling form of the qubit with its environment are. In general, ac- 
cording to the characteristic of the environment we set it a bosonic (or fermi) 
bath. The coupling of the qubit-bath can be supposed linear because the in- 
teraction between the qubit and the modes of the bath is very weak. Under 
these suppositions (linear coupling and short-time approximation) , V. Privman 
et al. 01 m investigated the decoherence of many kinds of open systems 
and we also investigated a concrete JCQ model in Ref.Pj. However, the non- 
linearly coupling forms of qubit-bath certainly exists 52 ■ In this Letter, we 
firstly introduce a model that JCQ nonlinearly couples with its bath and then 
use the short-time approximation to investigate the decoherence of the JCQ in 
the qubit-bath model. 



2 Nonlinearly coupling qubit-bath model 

Though the combined effects of the bath modes on the qubit may be large 
enough, the coupling of the qubit with each of the modes of the bath is very 
weak and the qubit-bath can be considered a linear coupling. However, due to 
the development of low temperature technique and fabrication technique of high 
quality quantum systems the nonlinearly coupling qubit-bath models can not 
be excluded in the devices of quantum. The linearly coupling qubit-bath model 
is analogous to the amplitude damping model. In this Letter, we introduce a 
nonlinearly coupling qubit-bath model which corresponds to the phase damping 
model We begin our investigation by considering a general open quantum 
system 

H = Hs + Hb+ Hi. (1) 

Here, Hs is the Hamiltonian of the investigated quantum system; Hb is the 
Hamiltonian of the bath; and Hj describes the interaction of the system and 
the bath. The bath is traditionally modeled by a large number of uncoupled 
harmonic oscillators (with ground energy shifted to zero), 

HB=J2''A^r (2) 

i 

If the interaction is linear, in the rotating-wave approximation the interaction 
Hamiltonian is 

if/ = A^.9,6t+At^5;6„ (3) 

j 3 

with the interaction constants gj . If the investigated system is a two-level system 
then A ^ (T-, ^ cr+. The short-time decoherence of this model has been 
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investigated in Ref."7". Suppose the bath be ohmic then 



(w) D {uj) = T]uj exp [uj/uJc) , (4) 

where D {uj) is the density of states of the bath and ojc is the cutoff frequency. 
For the model of JCQ couphng with the bath (see Fig.l of Ref . [T^ ) . the damping 
coefficient can be calculated with 



Rq \Cj 



where Rq = {2ef /h is the (superconducting) resistance quantum. When 
we investigate the short-time decoherence of this model, we can directly use 
{lo) D {oj) and do not need solving D (w) then g (u) . 

As well known that the interaction Haniiltonian of a oscillator in the phase 
damping bath is jTO] JT] 



which is a nonlinear coupling. If we replace the oscillator with a two-level 
system, namely, — + cr+, a — > cr_, then the total Hamiltonian becomes 



where g'^ = gj/'i, and oj'^ = g-j -I- ojj. This Hamiltonian models a kind of non- 
linear coupling of the bath to the exposed qubit. This bath can be considered 
resulting from the fluctuations of the Josephson nanocircuits and the surround- 
ing circuits. A more general model which includes other sources of decoherence, 
such as quasiparticle tunneling, fluctuating background charges and flux noise 
has been put forward ^31 dj- Suppose the bath is also ohmic then Eq.Q) will 
be held. But in this case, in order to investigate the decoherence of the qubit we 
must find out g (uij) through solving D (ujj). We suppose that the initial state 
of the environment state is G which is the product of the bath modes density 
matrices 



Trj (e-Z^A'^O ' 

where Mj = cojb^^bj. Because the degrees of freedom of the bath are much larger 
than that of the qubit, in the short time we can take the bath's states unchanged 
and consider the bath keeping in their thermal equilibrium state. We take the 
bath to be a mix of infinite harmonic oscillators at temperature T. Thus, the 
density of state of the bath reads 

exp('U) 
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for kT <^ hwj. Here (3 = l/kT, k is the Boltzmann constant. By using Eq.Q) 
we have 



5 K) = /^===exp(wj/2wc) 
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~ ^%/^exp(tjj/2a;c) • (10) 

3 Short-time dynamics of the qubit in the non- 
hnear qubit-bath model 

In this section we will investigate the general expression for the time evolution 
operator of the qubit in the nonlinear qubit-bath model within the short time. 
The definition of the short-time concept has been stated in Ref. [H] |7] The 
evolution operator is 

In the following we set t = t /h. Due to non-conservation of Hg in this system, the 
evolution operator cannot be in a general way expressed as ^-^Hst^-i{Hi+HB)t _ 
But in the sort-time approximation, the operator can be approximately ex- 
pressed as cni 

It has been proved that the expression is accurate enough for the time being 
short to the characteristic time. So the density matrix elements of the reduced 
density matrix p (t) in the basis of operator Hs can be expressed as 
p,„„ = Trs (^„Je-*^=*/2e-*(^^+^«)*e-'^=*/2i?(0)e^^=*/2e'(^^+^«)*e'^=*/2|^^j^ 

(13) 

Here, we suppose that the initial state of the system is R{0) — p (0) (8) 6, where 
p(0) is the initial state of the qubit and '9 — Yi^j- using the completeness 
relation of the eigenstates of Hg, i'l — have 

i=o,i 

where Xj [j = 0, 1) and \ipj) are eigenvalues and eigenstates of Hamiltonian Hg. 
Similarly, we have 

e^^^"'^"-^' = ^E^exp 1^ [±^ (g^X, +^) h]h,t] | \^^) (^^.| , (15) 

where Xj U — 0^ 1) ^-nd (j = 0, 1) , are the eigenvalues and eigenstates of 
operator cr^. So we have 

a=0,l 
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(3=0,1 p,g=0,l 
xPpjO)ne, E e**'" 1^^) (^J 



i/=0,l 



(16) 



namely, 

Pmn (t) 



a,/3,^,';,p,g,/i,^'— 0, 1 

{<Pa\ V'i) (V-il (P/s) (V5/3| 'Pp) {fql '^u) {<fi,\ (V'.l f^) 



{"Pj <^n) Ppg (0) 



(17) 



Taking the j term in Ttb [...], and denoting it by TVb [...] •, we have 



Trs [...]^. = TvB {exp [-i (^^Xj + ^) exp [i {g]x, + w^) } 

= Trs jexp -i [g'^x^ + Uj) b]bjt exp (^-f3ujb]bj^ 

[Tt exp (-/3a;,-6]6,)] ' exp [i {g]x, + ^) b]b,t\ } . (18) 



Because 



exp 



xbjbjt 



exp 



(e--l)6>, 



(19) 



where : : denotes the normally ordering form. By using the relationship of the 
completeness of the coherent states 



? = 



(20) 
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we have 
Ttb [. 



± exp \^(^e-<s>i+^'^y - l) b]bj 



exp 



-'-'-l)b]b, 



exp 



(e*(»>e+"^)* - l) b]bj'\ : |z) 



1 fd'z fd'i rd\ ( 2 ,„2 , , 

where Zj = (l — e^'^'^j ) ^ . By using the formula 

/d"^ z / 
exp ya \z\ +bz + cz* + fz^ + gz 

-abc — b'^g — c^f 



: exp 



we can obtain the resuh of (t) as 

Pmn (^) ~ 



(21) 



(22) 



^ g«i(A^ + A,^-Aa-A;3) 

a,/3,^,';,p,g,/i,i/— 0,1 

('/'ml <Pa) (^al V'^) (V'el ^ Vp) Wq\ V'^) {f^ V',-) 

It should be noted that it is not necessary to take N too large in the following 
numerical simulations because the ujj increase with j and the e~^'^^ decrease 
with ujj exponentially. By using Eq. ]2'6\i we can investigate various characteris- 
tics of the nonlinear qubit-bath systems. By using it, in the following we shall 
study the short-time decoherence of the JCQ. 



4 Measure of decoherence for a qubit 

In order to investigate the decoherence of the qubits, one must designate the 
measure to be used. There are many choices for this purpose. When the evo- 
lution time is very long, the qubit interacting with the large bath will fall into 
the thermal equilibrium at temperature T. In this case, Markovian type ap- 
proximation can be used to quantify the decoherent processes and it usually 
yields the exponential decay of the density matrix elements in the energy basis 
of the Hamiltonian Hg. In this time scale the measures of entropy and the first 
entropy can be used to quantify the decoherence. But the decoherence of the 
qubit gates cannot be characterized by this methods because the time of the 
elementary quantum gate operations are much shorter than the thermal relax- 
ation time. It has been shown that the norm ||(t||^ is useful for describing the 
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decoherence of the short-time evolutions. In this section, we will review the 
measure according to We set a is the deviation operator defined as 



a{t)=p{t)-p^it), (24) 

where p (t) and (t) are density matrixes of the "real" evolution (with interac- 
tion) and the "ideal" one (without interaction) of the investigated system. ||cr||_>^ 
is defined as 



A = sup ( ;.' I " ) • (25) 



For a qubit, the norm can be given by 



CTiol' + kiil'. (26) 

In general, the norm increases with time, reflecting the decoherence of the 
system. However, it is oscillating at the system's internal frequency. Thus, the 
decohering effect of the bath is better quantified by the maximal operator norm, 
D (t) which is defined as 

D{t)^snp{\\a{t,p{0))\\^). (27) 

P(0) 

It has been shown that the measure can rightly describe the short-time deco- 
herence of qubits. In the following we will use this measure to investigate the 
short-time decoherence of JCQ in nonlinearly coupling JCQ-bath model. 



5 Short-time decoherence of the JCQ 

In this section we firstly review the JCQ model. The single JCQ Hamiltonian 
is [IH] 

Hr ^ Ech{n~ngf - Ej cos ip. (28) 

Here, Ech — 2e^/ (Cg + Cj) is the charging energy; Ej — ^t^I^jlix is the Joseph- 
son coupling energy, where /c is the critical current of the Josephson junction, 
$0 — hc/2e denotes the flux quantum; Ug — CgVgjle is the dimensionless gate 
charge, where Cg is the gate capacitance, Vg the controllable gate voltage. When 
the Josephson coupling energy Ej is much smaller than the charging energy Ef-h , 
and both of them are much smaller than the superconducting energy gap A, the 
Hamiltonian Ea. (|28fl can be parameterized by the number of the Cooper pairs 
n on the island as ^7] JS| 

Hr ^ ^^Ech {n - Ugf \n) (n| 

n 

^]^Ej[\n){n + l\ + \n+\){n\]\. (29) 
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When the temperature T is low enough the system can be reduced to a two-state 
system (qubit) because all other charge states have much higher energy and can 
be neglected. So the Hamiltonian of the system can approximately reads 

Hs = -\b,<j, - ^B,a^. (30) 

Here, Bz — Ech (1 — 2?ig) and Ex = Ej. Ea. (|3()|l is similar to the ideal single 
qubit model 18 , but it can be modulated only by changing one parameter Bz- 
However, changing the parameter Bz (through switching the gate voltage) one 
can perform the desirable one-bit operations. If, for example, one chooses the 
idle state far to the left from the degeneracy point, the eigenstate lose to |0) and 
|1) . Then switching the system suddenly to the degeneracy point for a time t 
and back produces a rotation in spin space jl8| . 

TT f'^J* \ ( cos^ isin^fi \ 

Uj = exp [-^ax j = ^^.^^ ^^^^ j . (31) 

This can be obtained by modulating the gate voltage and making Bz — Ech 
(1 - 2ng) = 0. Thus, Hs = -^cFx, and Eq.Q) becomes 

H = -^ax + azY^ g',b]b, + ^ Lo'^b]b,. (32) 
j j 

In the following, we will calculate the Her (t)||^ and the decoherence D (t) of 
JCQ in the system In the calculation, three pure initial states are chosen, 
they are |(/))q = (1, 0)"^ , (corresponding to points in the following Figs.); = 

("^' l) ' (corresponding to above lines); \4>)2 = (^yh^ V^) ' ('^^'^'^^^P^'^'^^'^S 
to below lines), [because it has been shown that evaluation of the supremum 
over the initial density operators in order to find D (t) , see Ea. (l27|l one can do 
over only pure-state density operators '6 ]. Here, the eigenstates and eigenvalues 
of Hs are 

M = ^(|0)-|l)), Ao-^, 

= ^(|0) + |1)), Ai = -|^. (33) 



The eigenstates and eigenvalues of az are 



iV-o) = |0>, Xo = l, 

iV-i) = |1>, Xo = -l- (34) 



Thus, by using Ea.H23|l we can obtain 



an = \ (Poo - Pii) (2 - W^i - W^) + \p^^ sin ^ [W^ - W{) , (35) 
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o-io = ^(Poo-Pii)e ^ (14^1-1^2) 

+ \p,, [2 (1 + e^*^-0 + (14^1 + W2) (e**^'' - l)] , (36) 

where 

= n ZTTT^^'^^ = n (37) 

j=Na j=No 

In the following, we numerically investigate the decoherence of the JCQ in this 
model. We choose Ej = 51.8/iei; according to jj^], and T = 30mK. Suppose the 
frequencies of the harmonic oscillators made up of bath are range from toiow — 1 
MHz to ojhigh — 20 MHz = ojc (corresponding A'o = 1, iV = 20). Thus, we can 
plot the norms ||cr||;y versus time t in Fig.l. It is shown that when the initial state 
is p{0) = |(/))qq {(I>\ , ||cr(i)||;^ becomes maximum and it equals to D (t) (plotted 
by points in the Figs.). We denote the low decoherence {D {t) < 10~^) time 
^id_ Pj.q^j Fig.lb we obtain t'/ ~ 8 x 10"^ x 6.582 x 10"i°s = 5.266 x lO'V- 
We denote the elementary gate operation time, the characteristic time . In 
this case, = h/Ej = 1.27 x lO^ps. Thus, t'/ < t^. It shows that at the 
temperature T = SOmK the present setup of the JCQ cannot be taken as the 
qubit for quantum computation because it do not satisfy the DiVincenzo low 
decoherence criterion [Q. However, if we make the temperature T be lowered, 
the decoherence will decreases. For example, setting T — Q.3mK, and keeping 
Ej = 51.8iiev, we can plot the norms versus time t in Fig. 2. In this case 

we can obtain tl^ Ki2x 10"^ ^ 6.582 x 10-^°s = 13.16ps and tl=t{^ I2.7ps, 
where tif > tf. It is shown that when the temperature decreases to T = 0.3mK, 
within the whole time of elementary gate operation, D (t) < 10~*. Theoretically, 
in the lower temperature the JCQ in nonlinearly coupling JCQ-bath becomes 
an ideal qubit for making quantum computer because it satisfy the DiVincenzo 
low decoherence criterion. 



Fig.la, Fig.lb 



Fig.l: Norms W^W^ versus time t, (a) in a longer time comparing to the elemen- 
tary gate operation time; (b) in the low decoherence (D [t) < time. Here, 
the points and lines correspond to different initial states (see the explanation in 
text), Ej — bl.Sfiev, T — 30mk and Nq — 1, N = 20. The unit of the time is 
6.582 X 10"i°s. 
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Fig.2a, Fig.2b 



Fig. 2: Norms \\o'\\^ versus time t, (a) in a longer time comparing to the ele- 
mentary gate operation time; (h) in the low decoherence (D(t) < 10~*j time. 
Here, the values and meanings of the parameters are same as Fig. 1 expect for 
T = 0.3mk. 

6 Conclusions 

In this Letter we firstly constructed a qubit-bath model where we supposed 
that the qubit nonlinear ly couples with its environment. Then using the model 
we investigated the short-time decoherence of the JCQ. The decoherence is 
described by the norm of the deviation density operator, || cr (t)!];,^ . We define 
the biggest ||cr(t)||;^ the quantity of decoherence D{t). In the previous paper 
P] by using the same setup of the JCQ and bath's parameters we calculated 
the decoherence of the JCQ, where the JCQ is supposed linearly coupling with 
its bath. In the linearly coupling model, the decoherence is not larger than the 
DiVincenzo low decoherence criterion so the JCQ can be taken as the qubit for 
quantum computations. However, the results of this Letter show that at the 
same experimental temperature T = 30mA:, when the JCQ nonlinear ly couples 
with its bath, it is not a ideal block to build the quantum computer because 
the decoherence is larger than the DiVincenzo low decoherence criterion. But 
when the temperature decreases to T = 0.3mA: the decoherence decreases to a 
endurable grade. It suggests that in order to make the JCQ become a block of 
quantum computing hardware developing the technique of low temperature may 
be important. It should be noted that this is a pilot study to the decoherent 
problem of JCQ. The frequencies of the harmonic oscillators of the bath modes 
are chosen tentatively. They may only be obtained by virtue of some ingenious 
experiments. A further numerical simulation shows that when the frequencies 
are larger than our setting the JCQ will adapt a higher temperature. But, no 
matter what they are, we can find out a proper temperature for the JCQ in low 
decoherence by using our method. 
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